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ABSTRACT 

Turbulent channel flow and homogeneous shear flow have served as basic building block 
flows for the testing and calibration of Reynolds stress models. In this paper, a direct 
theoretical connection is made between homogeneous shear flow in equilibrium and the log- 
layer of fully-developed turbulent channel flow. It is shown that if a second-order closure 
model is calibrated to yield good equilibrium values for homogeneous shear flow it will also 
yield good results for the log-layer of channel flow provided that the Rotta coefficient is 
not too far removed from one. Most of the commonly used second-order closure models 
introduce an ad hoc wall reflection term in order to mask deficient predictions for the log- 
layer of channel flow that arise either from an inaccurate calibration of homogeneous shear 
flow or from the use of a Rotta coefficient that is too large. Illustrative model calculations 
are presented to demonstrate this point which has important implications for turbulence 
modeling. 


‘Research was supported by the National Aeronautics and Space Administration under NASA Contract 
Nos. NAS1-18605 and NAS1-19480 while the second author was in residence at the Institute for Computer 
Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. INTRODUCTION 


Turbulence models have been calibrated and tested using a variety of benchmark tur- 
bulent flows among which homogeneous shear flow and channel flow have played a central 
role. Typically these two flows are treated as separate tests that are completely indepen- 
dent. Starting with the work of Launder, Reece and Rodi 1 , the pressure-strain correlation 
of turbulence - which forms a pivotal part of second-order closure models - has been cal- 
ibrated based on the equilibrium Reynolds stress anisotropies in homogeneous shear flow. 
An ad hoc wall reflection term is then added to the pressure-strain model to yield good 
predictions for the log-layer of fully-developed turbulent channel flow. There are several 
disturbing features about the resulting model: the wall reflection term plays an important 
role far into the interior of the channel and it depends in an empirical manner on the normal 
distance from the wall. The latter deficiency makes it virtually impossible to systematically 
apply second-order closure models to turbule nt wall-bounded flows in complex geometries 
containing sharp corners. This, as well as other near-wall problems, has impeded progress 
in the application of second-order closures to the turbulent flows of technological interest. 

In this paper, it is shown that a second-order closure model will yield the same equilibrium 
Reynolds stress anisotropies in homogeneous shear flow and in the log-layer of channel flow 
if the slow pressure-strain correlation is represented by a Rotta 2 type of return-to-isotropy 
model with a coefficient of one. Since experiments 3,4 indicate that the Reynolds stress 
anisotropies for these two problems are close to one another, it follows that if a second-order 
closure model yields good equilibrium values for homogeneous shear flow it will also yield 
good results for the log-layer of channel flow provided that the Rotta coefficient is not too far 
removed from one. Illustrative calculations will be presented for four independent pressure- 
strain models — which include the models of Launder, Reece and Rodi 1 , Shih and Lumley 5 , 
Fu, Launder and Tselepidakis 6 , and Speziale, Sarkar and Gatski 7 - in order to demonstrate 
this point. Some rather surprising results are obtained concerning the performance of these 
models in channel flow. In addition, a crucial compatibility condition for the turbulent 
diffusion coefficient in the transport equation for the dissipation rate is elaborated on. The 
important implications that these results have for the development of improved second-order 
closure models are discussed in detail. 

2. THEORETICAL ANALYSIS 

We consider incompressible turbulent flows for which the Reynolds-averaged Navier- 
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Stokes equations take the form 


dui dui 

-X- + UJ — 

dt dxj 


dp 

dx{ 


+ uV 2 


Ui 


drjj 

dx.j 


(1) 


V • u = 0 


( 2 ) 


where U{ is the mean velocity, p is the mean kinematic pressure, r,j = is the Reynolds 
stress tensor, and v is the kinematic viscosity of the fluid. Here, the Einstein summation 
convention applies to repeated indices, an overbar represents an ensemble mean, and a prime 
represents a fluctuating quantity. The Reynolds stress tensor is a solution of the transport 
equation 8 
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at high Reynolds numbers where 
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Vij = dr^ u ' iU ' jU ' k + p ' u ' Sjk + p >u 'M) 

are, respectively, the pressure-strain correlation, turbulent dissipation rate, and turbulent 
diffusion term; Dj Dt = d/dt-\-u- V denotes the mean convective time rate and Kolmogorov’s 
assumption of local isotropy has been invoked. 

The two equilibrium turbulent shear flows to be considered are unidirectional with the 
mean velocity gradient tensor 


4 = stat » <4) 

where S = dujdy (see Figure 1). For homogeneous shear flow, S is a constant, whereas for 
the log-layer of turbulent channel flow, S — u T j ny where u T is the friction velocity and k is 
the von Karman constant (in more familiar terms, u + = (1 / n)ln ?/ + + 5 in the log-layer where 
u + = u/u T and y + = yu T ju). In channel flow, the mean convective terms are identically 
zero and within the log-layer, turbulence production equals dissipation (V = e) and, hence, 
the molecular and turbulent diffusion terms in (3) can be neglected 9 . Consequently, the 
anisotropy tensor bij = (r,j — | K )/2K and shear parameter S K/e (where K — is 
the turbulent kinetic energy) achieve constant equilibrium values in the log-layer that are 
independent of the boundary conditions. In homogeneous shear flow the molecular and 
turbulent diffusion terms in (3) are identically zero and each component of the Reynolds 
stress tensor grows exponentially at the same rate so that the anisotropy tensor b tJ and shear 
parameter SK/e achieve equilibrium values that are independent of the initial conditions 10 . 
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It is thus clear that the structural equilibrium in homogeneous shear flow and the 
production-equals-dissipation equilibrium in the log-layer of turbulent channel flow are each 
characterized by the constraints SK/s = constant and = constant. The latter constraint 
is equivalent to Dbij/Dt = 0, or 


~Dt 


(V-e) 


K 


( 5 ) 


where V = —t u S is the turbulence production. The substitution of (4) and (5) into (3), 
with vanishing turbulent diffusion terms, yields the equation 


K 



£ r *'2c T j 2 r , n,j 

SK = “ T 6 " + SK 


l-Z-Sa 

3 SK 3 


(6) 


which is valid for an equilibrium homogeneous shear flow and for the log-layer of channel 
flow. We will consider second-order closure models where 


TTjj — n 


(S) 


+ n& 


(R) 


(?) 


and the slow pressure-strain correlation is represented by a Rotta 2 type of return-to- 
isotropy model 

n !f = ~ C 'J ( r b “ \ K 6 v) ^ 

whereas the rapid pressure-strain correlation Tl\*p is modeled in the general form 


nif = (9) 

Here, both the Rotta coefficient C\ and the fourth-rank tensor Mijki can be functions of b{j 
(see the Appendix). 

If we make use of the fact that 


V T\2 SK 

£ K £ ' 


( 10 ) 


along with (8)-(9), it is straightforward to show that (6) can be written in the equivalent 
form 


T ij T\2 

K K 
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( 11 ) 


where = A4,ji2(&) is specified by the pressure-strain model chosen. Hence, since 

Tij/K is directly related to it then becomes clear that a closed set of nonlinear algebraic 
equations for the non-zero components of the anisotropy tensor (&n, 612, ^22 an d 633) are 
obtained once P ft is specified. Since P ft — 1 for the log-layer of channel flow and P ft ~ 1.8 
for an equilibrium homogeneous shear flow, it is clear that the same equilibrium values will 
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be obtained for these respective problems only when the Rotta coefficient Cj - 1 (the limit 
in which the dependence of on R fe is eliminated in (11)). It is also clear that this result 
carries over to the more general tensorially quadratic return models of the form 7 

nff = -2 Cxebij + 6(C X - l)e (b ik b kj - U kl b kl 8 t ^j (12) 

where the coefficient C\ can be a function of the second and third invariants of This leads 
us to the central result of this paper: A second-order closure model will yield approximately 
the same equilibrium values for b t j in homogeneous shear flow and in the log-layer of channel 
flow provided that Rotta coefficient is sufficiently close to one. In the next section, model 
calculations will be presented to illustrate that with a Rotta constant C\ as large as 1.7 it is 
possible to obtain good results for both channel flow and homogeneous shear flow without 
an ad hoc wall reflection term. 

3. ILLUSTRATIVE MODEL CALCULATIONS 

Calculations will now be presented for four pressure-strain models: the Launder, Reece 
and Rodi (LRR) model 1 , the Shih-Lumley (SL) model 5 , the Fu, Launder and Tselepidakis 
(FLT) model 6 , and the Speziale, Sarkar and Gatski (SSG) model 7 (see the Appendix for more 
details on the models). The equilibrium values corresponding to these models are obtained 
by substituting a given pressure-strain model into (6) and solving the resulting nonlinear 
algebraic equations numerically after (10) is made use of to eliminate SK/e. For channel 
flow, Vfe is set equal to 1 whereas for homogeneous shear flow, Vje is taken to be 1.8. In 
Table 1, the equilibrium Reynolds stress anisotropies and shear parameter SK/e obtained 
from the various models are compared with the experimental data of Tavoularis and Karnik 3 
for homogeneous shear flow. Several observations concerning these results are noteworthy: 
(a) the SSG and FLT models are, by far, in the best agreement with the experimental data 
for homogeneous shear flow, (b) the LRR model does not do as well since it was calibrated 
based on the older and less complete experimental data of Champagne, Harris and Corrsin 11 , 
and (c) the SL model performs the worst since, in its calibration, homogeneous shear flow 
was not directly accounted for. In Table 2, the corresponding model predictions for the 
log-layer of channel flow are compared with experimental data 4 (here, an average is taken of 
the log-layer values which vary somewhat with y + ). Apparently, only the SSG model yields 
equilibrium values that are in close range of the experimental data. The FLT model - which 
performs well in homogeneous shear flow - does not do quite as well in channel flow. This 
is a direct consequence of the theoretical result derived in the previous section. If a model 
yields accurate results in homogeneous shear flow, good results will automatically follow for 
the log-layer of channel flow provided that the Rotta coefficient is sufficiently close to one. 
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In the SSG model, the Rotta coefficient C, = 1.7 is sufficiently close to one to guarantee 


(Ci-l)IIMHMl (13) 

mTHv/e 

for all i and j where || • || is any suitable norm (this is a sufficient condition, that follows 
directly from (11), which guarantees that results for kj in homogeneous shear flow and 
channel flow will be close to one another as indicated by experiments). On the other hand, 
due to its nonlinear dependence on the invariants of the Rotta coefficient C\ « 3 for 
the FLT model which explains why the normal Reynolds stress anisotropies in channel flow 
differ by as much as 25% from their counterparts for homogeneous shear flow. The same is 
true for the SL model since its Rotta coefficient C, is approximately 5 in homogeneous shear 
flow (however, unlike the FLT model, the SL model renders inaccurate predictions for both 
homogeneous shear flow and channel flow). The LRR model has a sufficiently small Rotta 
coefficient C\ « 1.5 so that the deviations between its predictions for b tj in homogeneous 
shear flow and in channel flow are not fatal. The problem with this model is that it was 
not optimally calibrated for homogeneous shear flow - a deficiency that is tied to the fact 
that this model was developed before the more accurate experimental data became available 
which clearly indicated that production exceeds dissipation. In the calibration of the LRR 
model, the production was set equal to the dissipation for homogeneous shear flow 1 . 

Some comments are in order concerning how these results compare with the more de- 
tailed model calculations of homogeneous shear flow by Speziale and co-workers 7 ’ 10 ’ 12 and 
the recent systematic calculations of channel flow by Demuren and Sarkar 13 . For these more 
complete calculations, the Reynolds stress transport equation (3) must be supplemented 
with a modeled transport equation for the turbulent dissipation rate e which is typically 
taken to be of the form 1 


De 

Dt 


r s v c e2 + — 

C ^K V ' Ct2 K + dx i 



(14) 


where C cU C e2 and C £ are constants whose values vary from one model to the next. For 
homogeneous shear flow, the diffusion terms in (14) vanish. Since, DK / Dt — V — £ for any 
homogeneous turbulence, it then follows that an equilibrium state is achieved where 

V _ C& - 1 (15) 

£ C. i-l 

in the limit as t - oo (see Speziale and Mac Giolla Mhuiris 10 ). Hence, the equilibrium values 
for the various models given on Table 1 are identical to those that would be obtained from full 
Reynolds stress transport calculations using the model (14) with ( C e2 - i - 1) = L8 ’ 

Since most of the models do not employ precisely the same values for C c \ and C e2 , there 
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are some small differences between the equilibrium values displayed in Table 1 and those 
published previously 7 ’ 10,12 . However, the calculations presented herein for homogeneous shear 
flow actually form a more objective basis for the comparison of Reynolds stress models since 
V / £ \s set to a common experimental equilibrium value and the calculations are then freed 
from dependence on the model chosen for the turbulent dissipation rate. 

There is also a compatibility relation for the log-layer of channel flow that needs to be 
discussed. Since in the log-layer du + (dy + — e + = I//cy + and as well as K are constant, 
it follows that 

8(C €l - C c2 )b\ 2 


a = 


K 2 (2b 2 2 + |) 


(16) 


for the modeled dissipation rate equation to be consistent. Full Reynolds stress calculations 
of channel flow with models that satisfy the consistency constraint (16) will be in close 
approximate agreement with our calculations. The minor differences between the equilibrium 
values given in Table 2 based on our log-layer analysis and those obtained by Demuren and 
Sarkar 13 based on full Reynolds stress calculations are due to turbulent diffusion effects and 
the fact that some of the models considered herein violate constraint (16). Since C c2 — C e \ 
is in the range of 0.40 - 0.45 for most of the commonly used models, it follows that in 
order to yield a von Karman constant of « = 0.41 (with the approximate log-layer values 
of b\ 2 —0.15 and b 22 « —0.14), the value of C c chosen should be in the range of 0.16 

- 0.18. This constraint should be made use of more carefully in the future formulation of 
second-order closure models. 

Finally, some comments are in order concerning the wall reflection term that is added 
to many pressure-strain models in second-order closures to yield acceptable predictions for 
the log-layer of turbulent channel flow. Typically, the wall reflection correction is of the 
general form 1 

n .” = [c»,f (*» - §**,) + ^ (17) 

where IT^- ^ is directly related to the rapid pressure-strain model in the absence of walls, y is 
the distance normal to the wall, and C w i and C w2 are empirical constants. Since 


A " 3/2 r , +3/2 

= kK +3/2 k2.5 (18) 

ey v ' 

in the log-layer, and since is typically chosen to be in the range of 0.1 - 1.0, it follows 
that the wall reflection term makes a significant contribution to the slow pressure-strain 
correlation (this needs to be the case for many pressure-strain models due to their poor 
performance in channel flow as shown in Table 2). The problem with this is clear. At high 
Reynolds numbers the log-layer extends far into the interior of the channel. To have an 
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ad hoc correction - that depends on the normal distance from the wall - play a significant 
role far into the interior of the fluid is dangerous. It seriously diminishes the possibility of 
applying these models in complex geometries with corners where the normal distance y from 
the wall is not uniquely defined. 

4. CONCLUSIONS 

A direct theoretical connection between the log-layer of turbulent channel flow and ho- 
mogeneous shear flow in equilibrium has been established. These flows have traditionally 
been treated as being independent tests since in the former flow there is a production-equals- 
dissipation equilibrium, with bounded turbulent kinetic energy and dissipation, whereas in 
the latter flow, production exceeds dissipation so that the turbulent kinetic energy and dis- 
sipation rate grow exponentially with time. However, both flows have a common theoretical 
thread that connects them: the anisotropy tensor 6,y and shear parameter SI\/e achieve 
equilibrium values that are independent of the initial/boundary conditions. It was shown 
that in the limit as the Rotta coefficient goes to one, a second-order closure model will yield 
the same equilibrium values for 6,y in the log-layer of channel flow and in homogeneous shear 
flow. Furthermore, it was demonstrated that with a Rotta coefficient (7i as large as 1.7 
- which is a value that allows for the collapse of a significant range of return to isotropy 
data 7 - a model that was calibrated to yield good equilibrium values for homogeneous shear 
flow (the SSG model) also performs well in the log-layer of channel flow without ad hoc 
corrections. Hence, it appears that a model can be calibrated to perform well in both flows 
provided that the Rotta coefficient is not too far removed from one. 

The results obtained in this study have important implications for turbulence modeling. 
It is rather disquieting how poorly many of the currently popular second-order closure mod- 
els perform in the log-layer of turbulent channel flow. These deficiencies have their origin in 
two major sources: an inaccurate calibration of the model for homogeneous shear flow or the 
use of a Rotta coefficient that is too far removed from one (a state of affairs that has arisen 
from the introduction of an empirical nonlinear dependence of C\ on the invariants of 
The introduction of an ad hoc wall reflection term to alleviate this problem has seriously 
inhibited the ability to apply second-order closure models to turbulent flows in complex 
geometries. Since turbulent channel flow is dynamically similar to a two-dimensional equi- 
librium turbulent boundary layer - which forms a cornerstone for many practical engineering 
applications - it is crucial to get this flow right without ad hoc corrections that make the 
model geometry-dependent. The results of this study clearly show that it is possible to do 
this. More attention needs to be paid to this issue in the future if second-order closure 
models are to have an impact on the calculation of complex wall-bounded turbulent flows. 
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APPENDIX 


The detailed form of the pressure-strain models considered in this paper are as follows: 
Launder, Reece & Rodi Model 


n,j — —2C\eb{j + —A ~S x j -f C 2 K (bikSjk + b jk Sik 

-§6«3«*i) + C 3 K(b,kW jk + b :k w ik ) 

where 

7 ? _ 2 f @ u i @ U j\ TT7 duf\ 

ij ~ 2 \dxj + dxj ’ Wii ~ 2 [dxj dx,J 

C\ = 1.5, C 2 = 1.75, C 3 = 1.31 

Shih & Lumley Model 

n»j = —fiebij + -K'Sij + 12a 5 K {bikSjk + bjkSik 

— 2 bkiSkib + -(2 - las)K(b ik W ik + bjkWik) 

4 _ _ 

+^•(6,6 "h bjlbi m ~Sim 2>b{fcS kibij 

-3 bkiSkibij) + -K(biibi m Wj m + b 3l b lm W tm ) 

where 

0 = 2 + j .11 j \J~R^ t ){12 / \fRe t + 80.1 /n[l + 62.4(-// + 2.3 III)}} 


F = l +91 1 + 21 1 II 
77 2 b,j b t j , III — -bijbjkbki 





(A2) 

(A3) 


m 


(A5) 

(A6) 

(A7) 

(A8) 

(A9) 
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Fu, Launder & Tselepidakis Model 


— ybklbkifiij 

F—h’Sij + 1.2 A (^b%k$jk + bjkSjk ~ 

9fi 4 

+Y \K{b ik W jk + bjkWik) + -KfabuSj H 

~\~bj k b k i~Sn 2b ik S kl b,j — 3 b k i~$ k i b t j ) 

+\k (b{ k b k iW ji + b ih b u W«) - jK \siI(b ik W jk 

+ bjkWikj + 12 {bikbkiW Imbmj + bjkbkffi 

where 

ft = 120 IIF l/2 + 2 F 1/2 -2,/3 2 = 144 UF 1/2 


h|j — 0\sb{ 


+ ft£ {bikbkj 


(• 410 ) 


(411) 


Speziale, Sarkar & Gatski Model 


n tj = -{20,8 + c;v)b v + c 2£ (b ik b k3 - h kl b kl s^ + (c 3 - cj//J)* 3 y 

+C 4 K (bi k ~Sj k + bjk^ik — yb k [S k [ 8 t j^ + C 5 K(bi k W ) k + bj k Wi k ) 


where 


Cl = 1.7, c;- 1 . 80 , C 2 = 4.2 
c 3 = i c; = 1.30, C 4 = 1.25 

5 

C; = 0.40, Ih = b„b,, 


(412) 


(413) 

( 414 ) 

(415) 


11 


Equilibrium 

Values 

LRR 

Model 

SL 

Model 

FLT 

Model 

SSG 

Model 

Experimental 

Data 

b u 

0.152 

0.120 

0.196 

0.218 

0.21 

b 12 

-0.186 

-0.121 

-0.151 

-0.164 

-0.16 

622 

-0.119 

-0.122 

-0.136 

-0.145 

-0.14 

b 33 

-0.033 

0.002 

-0.060 

-0.073 

-0.07 

SK/e 

4.83 

j 

7.44 

5.95 

5.50 

5.0 


Table 1. Comparison of the model predictions for the equilibrium values in homogeneous 
shear flow (V je — 1.8) with the experimental data of Tavoularis and Karnik 3 . 
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Equilibrium 

Values 

LRR 

Model 

SL 

Model 

FLT 

Model 

SSG 

Model 

Experimental 

Data 

bn 

0.129 

0.079 

0.141 

0.201 

0.22 

b\ 2 

-0.178 

-0.116 

-0.162 

-0.160 

-0.16 

b 22 

-0.101 

-0.082 

-0.099 

-0.127 

-0.15 

b 33 

-0.028 

0.003 

-0.042 

-0.074 

-0.07 

SK/e 

2.80 

4.30 

3.09 

3.12 

3.1 


Table 2. Comparison of the model predictions for the equilibrium values in the log-layer of 
turbulent channel flow {Vje = 1) with the mean experimental data of Laufer 4 . 
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z 

(b) Log-Layer of Channel Flow 



Figure 1. Schematic of the equilibrium turbulent flows: (a) Homogeneous shear flow and 
(b) log-layer of channel flow. 
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of channel flow provided that the Rotta coefficient is not too far removed from 
one. Most of the commonly used second-order closure models introduce an ad hoc 
wall reflection term in order to mask deficient predictions for the log-layer of 
channel flow that arise either from an inaccurate calibration of homogeneous shear 
flow or from the use of a Rotta coefficient that is too large. Illustrative 
model calculations are presented to demonstrate this point which has important 
implications for turbulence modeling. 
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